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^ The crossover from orthogonal to the unitary universality classes in the distribution of the 

anomalously localized states (ALS) in two-dimensional disordered conductors is traced as a function 
O ' of magnetic field. We demonstrate that the microscopic origin of the crossover is the change in the 

, symmetry of the underlying disorder configurations, that are responsible for ALS. These disorder 

configurations are of weak magnitude (compared to the Fermi energy) and of small size (compared 
to the mean free path). We find their shape explicitly by means of the direct optimal fiuctuation 
method. 
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^ • PACS numbers: 72.15.Rn, 71.23.An, 73.20.Fz 

c : . . . 

Statistical properties of the wave functions, ij{r), in disordered conductors have been the subject of an intensive 
theoretical study during the last decade Considerable attention was devoted to the so-called anomalously 

localized states (ALS). These states constitute the large-jV'P "tail" of the wave function distribution, where this 
distribution deviates strongly from the prediction of the random-matrix theory. Analytical results obtained to date 
are summarized in the comprehensive review [Q . The bulk of these results are obtained using different versions [^|j^] 
of the nonlinear cr-model. The latest works on the subject report on numerical tests of the theoretical predictions 
[0-|lO|. Calculations of the ALS density, based on non-linear a model yield exponentially small tails with 

^ different numerical factors in the exponent, depending on the universality class. The drawback of the non-linear 
Q ■ cr-model approach is that it fails to pinpoint the actual disorder realizations responsible for ALS. This is because 
O I within the machinery of the non-linear cr-model the averaging over the disorder is performed at an early stage of the 
calculation. In an attempt to overcome this deficiency, an alternative approach, the direct optimal fluctuation method, 
I I was proposed in Ref. pd| | . In two dimensions the optimal disorder configuration obtained in Ref . pd[ | had a shape of a 
^ . compact potential well with a radius of the order of the Fermi wavelength, supplemented by a periodic modulation at 
' large distances. An important accomplishment of Ref. [ pd| is the observation that, within a numerical factor in the 
04 , exponent, the non-linear-cr-models-based results can be reproduced by considering a simplest configuration yielding 
the ALS, namely a circular (in two dimensions) region of space surrounded by a high and thick potential barrier. 
Then the anomalously large values of the wave function, ip, can be characterized by a single, intuitively transparent, 
parameter 



r-|^o«tl7IV'^nr, (1) 



where "in" and "out" refer to the inner and outer boundaries of the barrier. As a direct consequence of the compactness 
^ , of the fluctuations, the probablity of their formation is insensitive to weak magnetic fileds, in apparent contradiction 
' to the cr-model predictions. On the basis of this contradiction the authors of Ref. ||ll[] questioned the ability of the 
^ , non-linear cr- models to handle the ALS. 

O ■ The goal of the present paper is to demonstrate that the direct optimal fluctuation approach captures the dependence 
of the density of ALS on the universality class. Moreover, it allows to trace the crossover between the orthogonal and 
unitary classes with increasing magnetic field. We restrict our consideration to the two-dimensional case. 

The sensitivity to a weak magnetic field within the direct optimal fluctuation approach is restored in two steps. 
First, instead of circle-shape fluctuations of Ref. |]ll| we consider fluctuations of ring-shape, schematically depicted 
in Fig. 1. This step is in accord with the paper by Karpov |]l^ , who has demonstrated that in three dimensions 
fluctuations of a toroidal shape, providing a given value of T, are more probable than spherical fluctuations, considered 
in Ref. Ring-shape fluctuations, having much bigger area, are sensitive to the magnetic field. However, this class 
of fluctuations alone does not exhibit a crossover between the orthogonal and unitary classes. To reveal this crossover, 
as a second step, we note that, on a ring, the states with angular momenta m and —m are degenerate in the absence 
of magnetic flux. This suggests that the fluctuations lifting this degeneracy cause the increase of V{Ep,T), which 
is the probability to find a state with energy Ep and a given value of T. The underlying reason for this increase is 
the following. A weak coupling, k, between the states \m) and |— m) leads to decrease in the formation probability 
of the fluctuation. This decrease is quadratic in k. On the other hand, the |m), | — m) splitting due to this coupling 
is proportional to \k\. This splitting results in a "gain" in T, which is also proportional to Thus, for small 
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coupling, the "gain" overweighs the "penalty" for creatin g th e coupling. The above argument is completely analogous 
to the Jahn- Teller effect [Q. Kusmartsev and Rashba [Q have employed this argument to demonstrate that the 
fluctuations responsible for the tail states of a four-fold degenerate semiconductor valence band are of ellipsoidal 
rather than spherical shape. Summarizing, in a zero magnetic field V{Ep,T) is determined by the warped ring-shape 
fluctuations of the random potential. In a finite magnetic field the degeneracy |m), |— to) is lifted; the warp does not 
lead to increase of V anymore. Then the fluctuations contributing to V are close to rings. Thus, in the statistics 
of rare events, the crossover between the orthogonal and unitary classes can be viewed as switching from warped to 
almost perfect rings. This is illustrated in Fig. 1(a). 

The above discussion suggests the following analytical steps. Firstly, the fluctuation V{r) of the white-noise random 
potential is presented as a sum of the angular harmonics V{r) = X^m ^('°)^™'^' where p = \r\. Secondly, in the 
probability 'P{V{r)} of the fluctuation V{r) determined by 



\^^'P\='^jdrV\r)=^^Yl ldpp\Vm{p)\\ (2) 
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only the angular harmonics Vo{p) and V±2mip) are retained. This is because the |to), | — to) coupling is provided by 
the harmonics ±2to of the random potential. In Eq. (H) t and v stand for the scattering time and the density of states, 
respectively. Then the system of equations relating the components Xm a-nd X-m of the wave function p^/'^ip{p^ (f)) 
takes the form 

= + {""^T^P) + y--^) X = E.X. (3) 

where M is the electron mass, x = (Xm , X-m) is a two-component wave function, the matrix V2 = diag{V2m, V-2m) 
is diagonal, and ai are the Pauli matrices. The effective potential in is a sum of a centrifugal potential and Vb(p) 



The turning point, pc, is determined by V^J^\pc) = Ep [see Fig. 1(b)]. As we will see below, the characteristic width, 
w, of the fluctuation Vb(p) and the barrier thickness, d, [see Fig. 1(b)] are related as w ^ <C Pc- Similarly to 
Ref. [Q, this separation of spatial scales allows a number of crucial simplifications: (a) since w ^ d, it follows that, 
within the barrier, yj^^(p) is dominated by the centrifugal term; (b) since d ^ pc, this term can be linearized within 
the barrier region as 

V^}''^ip) = Ep + eo^ , (5) 

where eo ~ h^m^d/AIp^ is the barrier height. Using the above simplifications the ratio T = \tpEiPc — d,)\'^ /l^pEiPc)]^ 
can be readily calculated 



It is seen from Eq. (||) that the condition pc ^ d can be rewritten as to ^ | lnT|. Equation (||) relates T to the 
energy level position Eq in the potential well centered around pc~d^w/2 [see Fig. 1(b)]. Thus the problem of finding 
the distribution of T reduces to the conventional problem p5|Jl^ of finding the most probable configuration of random 
potential which confines the energy level of a depth —Eq. In our case —Eq is the eigenvalue of the matrix Hamiltonian 
H, defined as 

^--i^^ + yoi^) + V.a., (7) 

where x = p ~ pc + d + w/2. The form of the Hamiltonian Eq. (0) follows from the system Eq. (^) and the condition 
w d which allows to neglect the change of the centrifugal potential within the potential well [see Fig. 1(b)]. 

The straightforward generalization of the approach of Ref. [ p^ , p^ to the matrix Hamiltonian (|^) implies minimizing 
the functional F{Vo, V±2m} = \ InT'l — X{x^Hx)- With an appropriate choice of the Lagrange multiplier. A, we obtain 
the following expressions for Vq and V±2m in terms of two-component wave function, x, 
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= X+X , V±2^ = X+^±X • (8) 

At this point it is convenient to introduce dimensionless variables as z = x{2M£qY/'^ /h and X = ^^^X- With new 
notations, the matrix Schrodinger equation, Hx = ^£oX^ reduces to 

X" + \ (x+x) X~\ (x+f^^x) -^.X - X - , (9) 
and probabihty (Eq. (H)) of the fluctuation, providing a given value of T, takes the form 

\\nV,n\ =7T^'^ Jdx {\Vo{x)\^ + 2 \V2„,{x)f) = {3/m)7TkFl\ InTIC™ , (10) 

where kp — {2MEfY^'^ /fi is the Fermi momentum, I — fikpT/M is the mean free path; dimensionless factor Cm in 
Eq. (O) is defined as 



Cm = i I dz 



3(X+X)'-(X+<^.X)' . (11) 



Equations @-(^ express analytically the Jahn- Teller physics discussed in the introduction. Indeed, this system has 
a conventional instanton solution X-m = and Xm = 2^^^ cosh~^ z, corresponding to the absence of warp (V±2m = 0). 
This solution yields Cm — 16/3. On the other hand, the Jahn- Teller-type solution Xm — X-m = (2/3)"'^/^ cosh~^ z 
leads to a smaller value Cm — 32/9, which corresponds to exponentially higher probability Vm- Thus, the final result 
for orthogonal case reads 

ln7'™ = -(2/3)^5|lnr| , (12) 

where g = kpl is the dimensionless conductance. Remarkably, InVm does not depend on the value of the angular 
momentum m, provided that m 3> |lnT|. In fact, the condition m 3> |lnT| justifies all the assumptions made in 
deriving the result Eq. (|l^). To see this, we write explicitly the potentials Vq and V±2m' 

3cosh'[xi2MsoV/yn] ■ ^'^^ 

Using the relation Eq = Ep (3| lnT|/2m)^^'^, it is seen that under the condition m ^ \ lnT| we have Vo(x) <C Ep, i.e. 
the potential well is shallow. 

Other assumptions used in the calculation are less restrictive. Indeed, the typical width of the potential 
well w ~ ?i/(Meo)^/^ ^ kp^{m/\\nT\Y^'^ . On the other hand, for barrier width d we have from Eq. (H) 
d= (/9c/2)(3|lnr|/2m)2/3. Since = m/kp, we have d = mi/3(3| lnT|/2)2/3/2fcF, so that u;/d - |ln"^r| < 1. We 
have also assumed that the change of the centrifugal potential within the potential well is much smaller than Eq. This 
change can be estimated as SQw/d ~ eo/| lnT| <C eo- 

We now turn to the case of a finite magnetic field, B. Due to the azimuthal symmetry of the problem the system 
of equations Eq. (^ relating the radial functions Xrn,X-m can be easily modified to 

^"^+2M(r^+4r^''"^' ^''^ 

where Ib = {ch/ eB)^/'^ is the magnetic length. As can be seen from Eq. (^4|) there exists a following hierarchy 
of magnetic field strengths. When Ib '3> m^^'^/kp, then the "confining" term p^/l% can be neglected. The latter 
condition is equivalent to Rp 3> Pc, where Rp = kplg is the Larmour radius for an electron with energy Ep. On the 
other hand, it is obvious that the term -iifi^m/ (2M/^), lifting the |m), | — m) degeneracy, becomes important when it 
is of the order of eo- This yields Rp ^ pcilnT /m)~'^/^ . Since | lnT| ^ m, we conclude that there exists an interval of 
magnetic fields which affect the potential well but do not affect the barrier. Condition eg ~ ?i^m/(M/^) determines 
the characteristic magnetic field for orthogonal-unitary crossover: B^") = ($o/27r)fc|,(3|lnr|/2)2/3m-5/3, where $o 

is the flux quantum. The remaining task is to find the dependence Cm{5), where 6 — B/Bq"^\ which describes the 
change of Cm between "orthogonal", Cm(0) = 32/9, and "unitary". Cm = 16/3, values with increasing magnetic field, 
6. At finite 5 the system Eq. (^) takes the form 

3 1 

x" + o (x^x) X - o (x^o-zx) ^zX - X = <5(1 - ^z)x, (15) 
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whereas Eqs. (p^,(|Tl|) retain their form. At small S <^ 1 the dependence Cm(<5) can be found by solving Eq. ( |T5| ) 
perturbatively and substituting x{S) iiito Eq. (11). This program can be carried out analytically, yielding 



C™(5)=C™(0)(^l + y^. (16) 

Our important observation is that perturbative expression Eq. (|l^ ) remains valid with high accuracy up to (5 = 1/3, 
when the "unitary" solution X-m = 0, Xm = 2-'^/^cosh~^ z takes over. The crossover behaviour Cm{S) is illustrated 
in Fig. 2. In the same figure we show the "orthogonal" solution of Eq. ( |l5| ) at crossover point. 

Overall, the direct optimal fluctuation approach yields in two dimensions jlnT'l ~ gjlnTj rather than jlnT'l ^ 
gln^ T in Ref. |l^. This is because the truly optimal fluctuations are close to rings rather than to circles. The ring 
area, Sm ~ PcW, being much bigger than the area of a circle ~ kp'^ in is the origin of a high sensitivity of rings to 
the magnetic field. On the other hand, due to the same relation, Sm S> kp'^, the rings are much more "vulnerable" 
to the perturbations caused by the harmonics with small angular momenta, for which the centrifugal barrier is low. 
To demonstrate this we write down the expression for the decay rate of quasilocal state in a ring 

lmE = 2nJ2{HVo\^i)' + J2{m\Vn\^i)A6{E-E^)^- + -, (17) 

where (m|l^|/x) are the matrix elements of the harmonics Vn{p) [see Eq. (^]; \fi) is the state of the continuous spectrum 
in the absence of random potential (the corresponding energy is E^^). For simplicity in Eq. (|l^) we have neglected 
the warp. Above we were looking for the fluctuations with anomalously large T. Note, that within a numerical factor 
T — Ti/eoTi. It is seen from Eq. (|l|) that our consideration is justified only if T2 ti. Within the optimal fluctuation 
approach harmonics with small n, contributing to h/T2, do not affect the principal exponent. They enter at the stage 
of the prcfactor calculation [ pT[ . A typical value of T2 is the scattering time t, which is much shorter than ri. The 
relation T2 ^ ti is satisfied only for sparse configurations, in which the harmonics Vn are suppressed within the ring. 
Therefore, our result for Vm should be multiplied by the probability, P, to find such a configuration. To calculate 
this probability we consider the distribution function of time T2: 



P{t2) - 1 5 (—-( drdr'S{r,r')V{r)V{r')]) (18) 
\ V'^s J J lv{v) 



where the expression for S{r,r') follows from Eq. ( [17|) 



Sir, r') = '^p,^^{p)^*^{p') ^ ^,{r)rM) 5 [E - E,) . 



Averaging over random potential V{r) in Eq. (y^ can be carried out explicitly. The result is expressed in terms of 
eigenvalues, A„, of the integral operator with the kernel S{r,r'): J dr' S{r,r')(f)n{r') = A„0„(r). With the accuracy 
of a numerical factor the result can be obtained qualitatively. This is because P{t2) is determined by the first Nq 
eigenvalues, which are almost equal to each other: A„ ~ Xq for n < Nq. The value Nq can be estimated as the number 
of squares with a side kp^ within the area of a ring Sm, i.e. Nq ^ pcwkjp. On the other hand, there is an exact sum 
rule ^„ A„ = 2ttl' for the eigenvalues A„. Since A„ « Aq-ZVo, we can find Aq. This leads to the following result for 
the distribution Eq. ([l^) P{t2) oc {vt / \oT2)^° oc {t/t2)^°. Then the sought product Pm = T-'mP{T2 ~ ^i), describing 
the probability to find ALS with "preexponential" accuracy, can be presented as 

\nVm = InT'™ - C2A^o ln(T / g) = --TTg\ InT] - €2^7^ ln(r /g), (19) 

•J In T 

where C2 '--^ 1 is a numerical factor. For concreteness we choose in Eq. (|l^) the orthogonal expression for InVm- The 
above expression for Vm allows to find the optimal value of the angular momentum, m. Indeed, the prefactor falls off 
rapidly with m, while the main term to the first order is independent on m. We note that the m-dependent correction 
to the principal term in InV, which originates from the corrections to the barrier potential Eq. (^, is of the order of 
I \nV\d/pc, i.e., it has the form —cig(lnT/m)^/'^\ lnT|, where ci ~ 1 is another numerical coefhcient. The sign of the 

correction is determined by the fact that, due to deviation of V^J^j}{p) from the linear form Eq. (5), preserving the 
value of InT requires an increase in the binding energy, eg- Thus, the correction is negative, and increases with m. 
Probability Vm is maximal for m = ruopt = co{g\ lnT|)^/^, where cq = (ci/2c2)^/^. The resulting expression for Vm 
with the improved accuracy reads 
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2 / q r' 

llnP^I^-TTffllnTI (l + ^ 

It is seen from Eq. ( po|) that the above consideration is vahd within the domain 1 ^ InT <C g. For the optimal m 
the ring radius is Pc — CQkp^[g\ lnT|)^/^ <C I, whereas the ring width is u; = c^^^k'p^{g/\ lnT|)^/^. 

Note in conclusion, that the function V{T) studied in the present paper describes the distribution of the ALS widths. 
Therefore, among asymptotical distributions of various characteristics of ALS (see Ref. [||), 'P(T) corresponds to the 
distribution of the relaxation times. The random potential fluctuations that determine this distribution are shallow 
compared to Ep- Whether these fluctuations emerge within cr-model-based calculations remains to be seen. 
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FIG. 1. (a) Illustration of the crossover between orthogonal and unitary universality classes. Equipotential lines of the 
warped and ring-shaped fluctuations are shown schematically, (b) Effective potential for the angular harmonics m. Optimal 
Vb(p) at (5 = (zero magnetic fleld) and at critical 5 m 1/3 are plotted with solid and dashed lines, respectively. 
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FIG. 2. Normalized logarithm of the ALS density is shown versus the dimensionless magnetic field, 5. The inset shows the 
'orthogonal" optimal wave function (xm, X-m) for crossover value 5 « 1/3. 
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